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EXPERIMENTAL AND THEORETICAL PERFORMANCE OF COAXIAL PLASMA GUNS* 
by Charles J. Michels and Albert E. Johansen 
Lewis Research Center 

SUMMARY 

Experimentally determined calorimetric exhaust efficiencies and exhaust velocities 
are compared with predicted values from a snowplow model analysis. Four gun geome- 
tries (employing argon propellant) are examined analytically and experimentally over a 
range of initial masses and voltages. Noncrowbarred and crowbarred results are de- 
scribed. Experimentally stable operation of a gun has been limited to mass distributions 
for which a simplified analysis limits maximum efficiencies to about 40 percent. Experi- 
mentally measured efficiencies for optimum propellant distributions are approximately 
one-half those predicted from the simplified analysis (the maximum experimental effi- 
ciency was 22. 2 percent). Experimental efficiencies agree somewhat better with peak 
efficiencies predicted by the complete theory. Experimental velocities agree approxi- 
mately with both forms of the theory, but are generally somewhat lower. 


INTRODUCTION 

During 1964-65, experiments with coaxial plasma guns (all powered by conventional 
capacitor banks) were conducted at the Lewis Research Center. Measurements of kinetic 
efficiency based on calorimetric techniques and of magnetic front velocity of the plasma 
sheet in the gun were made for several gun geometries with wide ranges of injected mass 
and initial bank voltage. 

In order to interpret the experimental results, correlate the data, and determine 
the parameters important for efficient gun design, a theoretical model was developed and 
analyzed. The initial description of the theory for this model is given in reference 1. 


Part of this report incorporates, in condensed fashion, portions of a paper pre- 
sented at the Second Annual Meeting of the AIAA, given at San Francisco, California on 
July 26, 1965 (ref. 1). 



Like many previous theories (refs. 2 to 5), it is based on the so-called snowplow 
model of mass sweep-up by the current sheet but allows for a variable initial mass dis- 
tribution and includes, as a special case, the slug model (refs. 6 to 9). * For conven- 
ience, the mass distributions used to calculate the theoretical results presented herein 
are limited to members of a one-parameter family ranging from a constant distribution 
to that of a slug. In addition to variable initial mass distribution, this analysis includes 
the effect of losses associated with wall drag due to ion thermal diffusion, electron heat 
conduction, ohmic losses, ionization, and radiation. The assumed model leads to three 
coupled differential equations: one for the charge on the capacitor (Kirchoff's law), one 
for the displacement of the discharge (the momentum equation), and one for the plasma 
temperature (the plasma energy equation). The equations were integrated numerically 
with a high-speed digital computer. ' 

Theoretical results for various gun geometries are given for a simplified version of 
the analysis and also for the case in which all the mentioned losses are included. The 
simplified version includes only ohmic and shock heating losses. These results are pre- 
sented for comparison with the experimental data for both the noncrowbarred and crow- 
barred modes of gun operation. 

The experimental data allow only limited comparison with the theory. The total in- 
jected mass and the initial bank voltages were varied in a systematic way for each gun 
geometry investigated, while the other quantities were held constant. Some argon data 
previously reported by Michels and Ramins (ref. 15) have been included. The magnetic 
front velocities are measured in each gun annulus. No attempt is made to describe the 
plasma structure after it leaves the gun muzzle. This structure depends on the mode of 
operation of the gun and is described for a particular case by Michels (ref. 16). Ion ve- 
locity distributions in the exhaust, as described for particular cases in references 17 
to 20, were not examined in this work. 

An external crowbar system is described in appendix B by Fred F. Terdan. 

The meter -kilogram-second system of units is used throughout and the symbols are 
defined in appendix A. 


THEORY 

Three coupled nonlinear differential equations describe the capacitor-driven plasma 
gun system. The first two equations, the circuit and the momentum equations, involve 
functions of plasma temperature. The third equation is the plasma energy equation, 

1 Other models that have been analytically treated include the gas dynamic model 
(refs. 10 to 12) and the elastic piston model (ref. 13). The snowplow model has been 
justified to a limited extent by Hart (ref. 5) and Mather (ref. 14). 
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Figure 1. - Equivalent gun circuit. 


which is necessary to determine plasma temperature. The derivation of the theory pre- 
sented in reference 1 is reviewed in the following sections. 


Circuit Equations 


The equivalent gun circuit is shown in figure 1. Before crowbar occurs, 1^ is zero, 
and Ig = Iq = I. Kirchoff's equation for a circuit may be obtained by deriving an energy 
equation for the circuit, differentiating it with respect to time, and dividing it by the cur- 
rent. This approach is used in this analysis to account conveniently for the electrical 
energy losses associated with ionization and excitation of the plasma. The fraction of the 
ionization and excitation that is not produced directly by electrical energy is considered 
to be accomplished at the expense of plasma thermal energy. This fraction, denoted by 
P(T), is a function of the plasma temperature. It is assumed that the plasma is always 
100 percent ionized and a constant fraction g (assumed unity herein) of the ions are in an 
excited state corresponding to some single effective quantum level estimated from spec- 
troscopic observations and knowledge of excited states. 

With these assumptions, the electrical energy absorbed in ionizing and exciting par- 
ticles that are swept up by the discharge current sheet and in reionizing and reexciting 
particles that have diffused to the walls may be written as 



_d 

dt 



+ T(T) 


dt 


where + g$ e , M is the mass swept up by the current sheet, m. 


is the ion mass, 
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and r(T) is the diffusion rate of ions to the walls. The symbols and are, in 
turn, defined as the effective ionization and excitation potentials. The electrical energy 
lost by radiation (i. e. , by maintaining the excitation of ions in the discharge) may be 
written as 



where <2> 0 = g(<f> /t ) is the excitation loss rate parameter and t is the mean lifetime 
of the excited state. 

The remaining terms in the energy equation for the circuit are familiar and should 
require no explanation. The complete circuit energy equation is 


\ cv ? - \ cv2 < L P + l d + L ' x > l2 + C (R P + Rp)l2 


dt 


[ F *f dt+ f [ 1 - p<T) ]' 

Jo Jo 




— /— )+r(T) 

dt \rnj 




where 


x 


= - L 1 ! 2 


is the x-component of the Lorentz force integrated over the discharge volume. Differen- 
tiation of the energy equation and subsequent division by I (where I = -C(dV/dt) produces 
the following circuit equation: 

-V + (L + L D + L'x) — + — L T — I + (R + R p )I 
P u dt 2 dt P P 


- L' — I + - [l - — (— \ + T(T) + (— 

2 dt I L J I 1 dt \m i J J \m i 


^ 1*2 “0 


( 2 ) 


The terms in this equation correspond, respectively, to the following voltage drops: 
condenser (negative), inductance, increase in new magnetic field volume, resistance, 


4 



work done on the discharge, and finally, ionization and excitation (radiation). Note that 
the third and fifth terms are identical; that is, the voltage drop associated with the in- 
crease in magnetic field is the same as the voltage drop associated with the mechanical 
work done on the plasma. When the terms are collected, the circuit equation may be 
written as 


V = (R„ + Rp)I + — 
p v dt 



( 3 ) 


Equation (3) has been derived for a single moving discharge. In typical gun opera- 
tion, a secondary discharge occurs which may be either moving or stationary. If station- 
ary, it is called a crowbar discharge. The resulting circuit diagram is also illustrated 
in figure 1. Equation (3) is used prior to the time when crowbar occurs. When crowbar 
occurs, the crowbar "switch” may be assumed to close, and the following equations de- 
scribe the resulting circuit: 


V ~ R p X B + R P r G 


+ f t Vb + Vg + L ' xI G> 



o - + R c 1 c + dt ( L D X G + L C I C + L xI G^ 



Momentum Equation 

The snowplow model assumes that the discharge current sheet sweeps up the mass 
it encounters and carries it at sheet velocity. In addition, all the directed momentum of 
a particle striking the gun walls is assumed to be lost to the wall. The momentum equa- 
tion is therefore 


i( M f)i L,I °- m i r(T, f ,4, 
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which merely states that the time rate of change of plasma momentum is equal to the 
force on the discharge minus the momentum lost to the walls by ion thermal diffusion. 


Plasma Energy Equation 

The plasma energy equation is written for a frame of reference moving with the cur- 
rent sheet and is used to determine the static plasma temperature T that is required to 
evaluate T. In this frame of reference, each particle entering the discharge (either ini- 
tially or following contact with a wall) appears to be moving with a velocity -dx/dt and 
possesses a kinetic energy 1/2 m. (dx/dt) . The plasma energy equation states that the 
rate of thermal energy increase can be equated to the algebraic sum of the rate of heating 
due to thermalization of the kinetic energy of particles entering the discharge and the rate 
of cooling due to various assumed loss mechanisms. It can be written as 


d^ 

dt 


— kT(l + z) P 4 

2 K 


= Rplp + - — f—) 2 + - rn.T(T) f—Y - - kT(l + z)T(T) 
2 dt \dt/ 2 1 \dt/ 2 


- P(T)<$j 


A(M- ]+ r(T) 

dt mJ _ 



(-KVT)dS - - kT — 
2 e 


( 5 ) 


where the terms on the right side of the equation are respectively due to, ohmic heating, 
thermalization of the kinetic energy of the swept-up mass, rethermalization of the parti- 
cles rebounding from the walls, convective losses associated with particle diffusion, ion- 
ization and radiation losses, conductive losses of heat in the electron gas, and convective 
losses associated with the main electron conduction current. Wall drag effects are taken 
into account in the last term of equation (4) and the fourth term on the right side of equa- 
tion (5). Approximate expressions for T(T), f (-KVT)dS, and K used throughout 

... . •'wall 

this report are 

2kT 


r (X ) = I (m\ 1 /ikl 

2 \ m i/ (r o " r i> i m i 



(-KVT)dS = 7.77 ttK Ax 
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and 


k(kT) 5 / 2 


The expression for r is just one- fourth of the product of the thermal velocity, number 
density, and area, correct for an isotropic velocity distribution. The rate of heat loss 
by conduction in the electron gas was deduced from an analysis of a simplified situation, 
namely, steady- state conduction between isothermal slabs with uniform heat source den- 
sity. The thermal conductivity of the electron gas K is as described in reference 21. 



Nondimensional Equations 

The three nonlinear coupled differential equations (eqs. (3), (4), and (5)) completely 
describe the system consisting of capacitor, gun, and plasma. The fraction P(T) must 
also be specified and was assumed zero for all the calculations. In order to generalize 
the equations and facilitate their solution, the following nondimensional variables and pa- 
rameters are introduced: 


r = (L VC)" 1/2 t; Se= /-L-) L ; £ = jT 1 x; f = MM; y = Q' 1 Q 

\ L V M 0 U 


= 2k/iLC\ T . jf = i M • 6 = 


m.\ £ 



l3 r o ' MJ 


£; t/— 

Y L’£ 


V 1 


(, + g* ); „ = AM Jpa g* „ = 3 PH 
V m iV \mjtl t e 2 ev 0 


m,£ 


/ = 


= 2.85 



,2/ i 


tL*Cy 


i 4 ire. 


In A 


R 


( 6 ) 
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col ^ 


A mass density distribution in the x-direction M Q (df/dx) was so chosen that a family 
of distributions could be represented from a "slug" to a constant density profile. The 
function 


f = 1 - (1 - |) 1 /( 1_Q! ) 


(7) 


was used in which a varies from zero (constant density) to 1 (slug). 

Before crowbar occurs, the nondimensional system equations can be written as 
follows: 


Circuit equation: 

y + (^ p + &£))y + (^ p + # p )y + kr + + 



f*i + SfjS 1 / 2 ) + <p 2 f] = 0 


Momentum equation: 


Jf'k 2 +Jii =r 2 - 5JfkP 1/2 


( 8 ) 


(9) 


Plasma energy equation: 

(1 + z)J(t'kP + f&) = 2# p y 2 +±J( f*i + 6f/3 1/2 )i 2 - | (1 + z)6 Jtp 3/2 - f(3 7/2 - Sl\y |/3 


( 10 ) 


Not an explicit part of the system of equations but available from the calculations 
are the kinetic efficiency r\ 




- M(— 
2 \dty 


lev 2 

2 u 


= (±M 2 
\2 ) 


4=1 


x=l 


and the nondimensional velocity £ 


£ = _ Vgj /dx N 


(dr. 


l=i Ve \ dt 4=f 
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CO I TjH 


After crowbar occurs, there is a similar set of nondimens ional equations that are 
complicated as a result of the extra circuit involved. These are written as follows: 


Circuit equations: 

y+ jsy>+ £ B y - - (iJ + 4>) - Rp> - /' 1 |^^ 1 (f , i + 6f/3*/ 2 ) + 



( 11 ) 


(£/ + U) + (* D + se c )j + Rpi + R c / + if c y + R c y + ± [«P 1 (f’6 + 6f/3 1/2 ) + <p 2 i] = 0 

( 12 ) 


Momentum equation: 


M'k 2 + Mi = (-/) 2 - 6JTfi/3 1/2 


(13) 


Plasma energy equation: 


(1 + zU(f f |/3 + #) = 2Rp(-/) 2 + A u*f(f’i + Sf^ 1 / 2 )! 2 

- 1 (1 + z)&Jip 3/2 - //3 7/2 - S3 1 (-> ) | /3 (14) 


In equations (11), (12), (13), and (14), the additional complication of crowbar circuitry 
requires that extra nondimensional current variables be defined, namely, 


y = 


A B 


CV 0 /(L*fC) 1 / 2 


/ = 


CV Q / (L’f C) 1 ^ 2 


Theoretical Performance 

The nondimensional system equations were solved by numerical techniques for two 
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Figure 2. - Effect of plasma cooling on kinetic efficiency for geometry A. Non- 
crowbarred; nondimensional mass distribution parameter, 1/2. 


different cases, a simplified model and a loss model. Ohmic losses and losses due to 
shock heating are included in the simplified model. The loss model includes the plasma 
cooling terms accounted for in equations (8), (9), and (10). The effect of these cooling 
terms on efficiency is shown in figure 2 for a particular case (noncrowbarred, a = 1/2, 
geometry A, described on p. 20). The effect of wall drag on efficiency is dominant. If 
no plasma cooling mechanisms are considered, the resulting high temperatures of the 
plasma sheet enhance the diffusion to the walls, and the loss of momentum to the walls 
limits the maximum efficiency to 7. 1 percent. Some cooling of the plasma is necessary 
to increase the efficiency to near observed values, which generally are above the values 
shown by the lower curves of figure 2. Although such cooling is also an energy loss, it 
is not a momentum loss. In fact, it reduces the drag loss so that the overall effect is an 
increase in kinetic efficiency. Results for other configurations are presented in the sec- 
tion COMPARISON OF THEORY WITH EXPERIMENT along with the corresponding experi- 
mental data. 

The useful thrust of the device is the sum of the rate of momentum flux through the 
exit plane and the force equal to the product of exit pressure and exit area. The theoreti- 
cal calculation of efficiency has not taken into account the pressure effect. In the simpli- 
fied model, sample calculations incorporating the pressure term have shown that theoret- 
ical efficiencies increase a few percent. In the loss model, however, the plasma cooling 
has reduced the sheet pressure so that the effect is negligible. 

In reference 1 it is shown that the simplified model efficiencies were greatest for the 
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slug mass distribution, for Jl near 1, and near 0.1. Also, crowbarring (at either 

maximum current time or maximum energy time, that is, the time when the useful energy 
available for plasma acceleration after crowbar is a maximum) does not give a higher 
maximum efficiency than the corresponding noncrowbar case. Crowbarring at the time of 
maximum available energy was found to have higher efficiency for Jl greater than 1 . 0 
than either of the other two modes of operation. If for other reasons, such as stability or 
breakdown delay time, a gun were operated at Jl greater than 1. 0, crowbarring at maxi- 
mum energy would give the highest efficiency, other factors being held constant. 

Crowbarred loss model cases are not included in reference 1, but they have since 
been calculated, and the results are presented herein. Generally, the incorporation of 
losses in the crowbarred analysis shows that efficiencies at values above Jl equal to 1. 5 
are less than the corresponding noncrowbarred case. Previously, crowbarring was con- 
sidered advantageous for Jl greater than 1.0 (ref. 1). This advantage, however, is 
more than negated by the fact that losses are severe for Jl greater than 1. 0. 


APPARATUS 


Capacitor Bank 


The guns described in this report were energized by a capacitor bank consisting of 
11 identical sections switched simultaneously through cabling to the coaxial plasma gun. 
Each section consists of a 1. 1-microfarad capacitor, a GL-7703 ignitron switch, and 
three parallel coaxial cables. The 5000-Joule bank, which has a resistance of 3. 3 milli- 
ohms and an inductance of 14 nanohenries, is discussed in detail in reference 22. 
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Figure 3. - Coaxial plasma gun. 


To 

vacuum 

pump 


11 



Gun Geometries 


Figure 3 shows the coaxial gun and the manner in which it is connected to the bank 
cables. The electrodes of each of the guns are made of oxygen- free high- conductivity 
copper. At the breech of the gun is a Vycor insulator and a gun header assembly. The 
header connects the gun electrodes through a parallel-plate transmission line to the ca- 
pacitor bank cables. The header also served as a convenient place to connect an external 
crowbar switch (described in appendix B) during some of the experiments. 

A short- duration puff of argon gas was admitted into the evacuated gun annulus 
through ports in the center electrode by a fast-acting mechanical gas valve. The duration 
of the gas puff was set experimentally at 100 microseconds. The mass per puff was con- 
trolled by the gas pressure in the valve plenum. An air-operated hammer provided the 
controlled blow to actuate the valve stem. 

The dimensions and electrical characteristics of the four gun geometries investigated 
are given in table I. Geometries A and B used the same electrode radii (radius ratio of 
about 3), but the gas port location was changed. Geometry C employed gas ports at the 
breech, a glass baffle (to prevent gas reaching the breech to form a static secondary, 
i. e. , a crowbar discharge), a low inductance breech configuration, and a larger radius 
ratio (7. 5). Geometry D differed from geometry C only in its radius ratio (15) and 
shorter barrel. This design caused the gun current to be near zero when the plasma 
sheet was at the muzzle. In all experiments (except when geometry A was so operated 
that crowbar occurred at time of maximum current), a corona inhibitor ring covered a 
corner protruding in the breech. This delayed the time of crowbar. 


Instrumentation 

Bank charging voltage was measured with calibrated voltmeters mounted on the bank 
control console. Gun exhaust energy was measured as the temperature rise of a calo- 
rimeter cup located 30 centimeters downstream of the gun. The temperature rise was 
recorded on a millivolt strip- chart recorder. The calorimeter cup (15 cm diam.) was 
constructed to reduce ablation errors and to absorb the beam energy efficiently. 

Probes for the measurement of azimuthal magnetic field were inserted into the gun 
annulus, 11. 5 centimeters apart, staggered along the gun axis (four in geometries A, B, 
and C and two in geometry D), and immersed in the annulus about 0. 3 centimeter so as 
not to affect the plasma sheet. Time of flight of the magnetic front between two adjacent 
probe stations determined the front velocity. The magnetic probes had four turn coils, 
solenoidally wound about a 0. 050- inch- diameter coil form and encased in a Vycor protec- 
tor tube. The signals were transmitted via coaxial cables, terminated, passively inte- 


v 
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TABLE I. - GUN GEOMETRIES AND CHARACTERISTICS 


[Radius of gun center electrode, r^ radius of gun outer electrode, r Q ; effective gun length, i; gun length from initial discharge loca- 
tion to crowbar discharge location, parasitic resistance, R p ; crowbar resistance, R^,; parasitic inductance, L p ; crowbar in- 
ductance, L^; gun inductance from crowbar discharge location to initial discharge location, Lp. ] 




grated (10 /isec time constant), and recorded on an oscilloscope. The probes were cali- 
brated by immersion in a known transient magnetic field with dynamics similar to that of 
the experiment. 

Gun voltage, gun current, and the sequence of events were used to monitor the exper- 
iments. 


EXPERIMENTAL PROCEDURE 

The first step in the operating sequence for the experiments was to evacuate the gun 

— fi 

and test section to the low (10~ mm Hg) pressure region. A short duration puff of argon 
propellant was then injected between the electrodes. A controlled period of time was 
allowed for the gas to disperse between the electrodes and for the charged capacitor bank 
to be switched to the gun. After breakdown, the resulting plasma accelerated out of the 
gun and into the evacuated test section where it impinged on the calorimeter cup. In all 
reported data the center electrode was initially negative. 

For each geometry, a series of shots was made to determine the optimum time delay 
between insertion of the gas puff and application of the bank voltage. This optimum delay 
time (with respect to kinetic efficiency) was used in the remainder of the shots with the 
same gun. 

The following procedure was used to obtain each data point shown in figure 4: 

(1) The mass per puff was calibrated just before the series of shots. 

(2) The gun was run 10 times to outgas and condition the electrodes. 

(3) All data were gathered simultaneously and synchronously for five identical shots. 

(4) The mass per puff was recalibrated. 

The experimental data points on the figures are the average values of the five shots 
gathered. 

The variation of efficiency and magnetic front velocity with mass over the practical 
operating range of the device was determined for three different bank voltages: 15, 20, 
and 25 kilovolts. The experimental gun efficiency was determined by dividing the energy 
received in the calorimeter cup by the difference between the energy originally stored in 
the bank and the switch loss energy (ref. 22). The calorimeter cup cannot differentiate 
between the thermal energy and the directed kinetic energy in the impinging exhaust. The 
exhaust has been assumed herein to be primarily directed kinetic energy. This assump- 
tion is justified, in part, by the fact that the theoretical thermal energy component is 
small when losses are included in the analysis (as discussed in the section Theoretical 
Performance) . 
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(a) Geometry A, not crowbarred. (b) Geometry A, early internal crowbar. Theoretical model crowbarred at maximum current. 

Figure 4. - Theoretical and experimental performance of various plasma-gun geometries. 
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Figure 4. - Continued. 
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(e) Geometry B, noncrowbar red. 


(f) Geometry B, delayed internal crowbar. Theoretical model crowbarred at maximum energy. 
Figure 4. - Continued. 









COMPARISON OF THEORY WITH EXPERIMENT 


The experimental and theoretical energy efficiencies and velocities determined for 
each geometry are presented as a function of the mass loading parameter Jl in figure 4. 
Theoretical efficiencies are presented for both the simplified and loss models. 

Neutral gas-pressure profiles in the gun as a function of time were measured prior 
to the experiment for a geometry similar to geometry A (without the glass baffle) . From 
knowing the time delay for optimum experimental efficiency and the shape of the pressure 
profile at that time the mass distribution factor a was estimated to be approximately 
1/2. Pressure profiles against time for the other geometries were not determined but 
an a of 1/2 was used for calculating all theoretical efficiencies. 

The loss parameters evaluated for each geometry in determining the theoretical effi- 
ciency for the loss model are shown in table n. These parameters were mostly deter- 
mined from the configuration in table I and the general or estimated values indicated in 
table n. 

The experimental efficiency curves are plotted for a limited range of Jl values. 

For low values of mass inserted in the gun annulus (low values of Jl ), the statistical 
breakdown delay time (time between voltage application and ignition) increases. Estima- 
tion of Mq and a for the theoretical analysis therefore becomes difficult, and the data 
are not repeatable because of the nonuniform breakdown time. The experimental curves 
terminate when Jl is reduced to the point where breakdown delay is longer than 10 micro- 
seconds. For large values of inserted mass, the discharge becomes diffuse. The curves 

TABLE II. - LOSS PARAMETERS 

[Thickness of current sheet, Ax, 0. 01 m; effective ionization po- 
tential per ion, 4>., 1. 12x10" ^ J (70 eV); effective excitation 

1 _ -1 O 

potential per ion, $ e , 6. 42x10 J (40 eV); fraction of ions in 
excited state, g, 1. 0; charge number of ions, z, 2. 0; mean 
lifetime of excited state, t 0 , 2. 5xl0‘ 8 sec. ] 
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terminate when the steepness of the magnetic fronts measured in the gun annulus becomes 
so small that estimating the sheet velocity becomes difficult. 

The maximum theoretical kinetic energy efficiency for the simplified model (ref. 1), 
for a equal to 1/2, is 40 percent when Jl equals 1.2 and i? tot equals 0.092. Even 
though geometries C and D are most nearly like this model, a survey at this Jl value 
showed that none of the geometries had an efficiency greater than 22. 2 percent. 

Other practical factors not included in the analysis also have effects on the maximum 
kinetic energy efficiency. For example, for a gas propellant, the simplified model opti- 
mum solutions suggest insertion of the gas as a slug at the breech. However, very poor 
efficiencies were noted when this mass distribution was attempted experimentally (by a 
very short delay between injection and bank voltage application). Magnetic probes indi- 
cated that "spoking instabilities, ” which prevent the efficient sweep-up of gas, exist for 
these cases. These instabilities decrease the efficiency well below the theoretical value 
for a. near 1, while even the theoretical efficiency is poor when a is near zero. That 
optimum experimental efficiency occurs for intermediate values of a (e. g. , a = 1/2 for 
geometry A, internally crowbarred) is therefore understandable. The experimental effi- 
ciency for most geometries was significantly affected by initial bank voltage. This effect 
did not agree with the theoretical analysis, which predicted a much smaller voltage effect 
in the opposite direction. The small effect of bank voltage permitted the use of one bank 
voltage, 25 kilovolts, to determine the loss parameters of table n and the theoretical 
curves of figure 4. 

The experimental and theoretical efficiency and velocity are shown as a function of 
the mass loading parameter JL for noncrowbarred geometry A in figure 4(a) (p. 15). A 
Vycor baffle was placed just upstream of the gas ports in the gun annulus for this case to 
prevent gas from drifting toward the breech and eventually providing conditions for devel- 
oping a second static crowbar discharge. The mass distribution for this mode of opera- 
tion was not measured. Theoretical loss model efficiency curves are shown for a equal 
to 1/2 and 0 to emphasize the need for accurate knowledge of mass distribution before 
attempting precision comparisons of theory with experiment. 

Figure 4(b) (p. 15) is for the same geometry operated without the glass baffle so that 
it crowbarred internally. In this experiment, the breech outer electrode connection had 
a relatively sharp radius edge on which the crowbar discharge could seat. ' The gun crow- 
barred at close to maximum current time. The 20 kilovolt experimental data straddle 
the loss model curve for Jl greater than 2. 0 to within 2 percent. The 15 kilovolt data 
are in poorer agreement except near Jl equal to 3.0. 

The same geometry, operated crowbarred internally without the glass baffle but with 
a ring to inhibit crowbar breakdown on the sharp edge at the breech, provided the data 
shown in figure 4(c) (p. 16). The presence of the ring delayed the time of crowbar rela- 
tive to the previously mentioned case (closer to the case for crowbar at maximum 
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energy). The experimental efficiencies are within 3 percent of the loss model predictions 
for Jl greater than 0. 7. The efficiencies increased for the delayed crowbarring as com- 
pared with the other crowbar case and peaked at values of Jl closer to 1. Both these 
quantitative findings appear to be in agreement with the simplified and loss analyses. 

An external crowbar switch was connected across the header terminals of geome- 
try A. Figure 4(d) (p. 16) shows the data for external crowbar at maximum current time 
(experimentally determined as optimum time of crowbar). In order that the crowbar 
switch efficiently trap the energy in the gun, the crowbar switch impedance should be 
about one-tenth the bank system impedance. A crowbar switch with impedance lower than 
an already very low impedance bank system was physically impractical. External crow- 
bar was nevertheless attempted, because independent control of crowbar switching could 
be achieved. 

A glass baffle was placed just upstream of the gas ports to prevent internal crowbar- 
ring at the breech for the experiment. The theoretical loss-model efficiency curve is 
similar to the corresponding noncrowbarred curve of figure 4(a), because the external 
crowbar impedance is relatively large. The experimental efficiency data for the external 
crowbar (fig. 4(d)) are generally a few percent lower than the corresponding noncrow- 
barred data (fig. 4(a)). This difference is probably due to the external crowbar switch 
losses. The comments concerning mass distribution (discussed on p. 20) for figure 4(a) 
apply to figure 4(d) also. 

Figure 4(e) (p. 17) shows the data for geometry B operated noncrowbarred. A glass 
baffle was inserted upstream of the gas ports to prevent gas from drifting toward the 
breech and eventually to provide conditions for developing a crowbar discharge. Limited 
experimental data were available. 

In the experimental curves of geometry B operated internally crowbarred (fig. 4(f), 
p. 17), the voltage again causes some shifts in the efficiency curves. The crowbar oc- 
curred later than at maximum current time and is compared with theoretical values for 
crowbar at maximum energy time. This geometry operated at very low injected masses, 
and breakdown delay times beyond 10 microseconds made data gathering impractical for 
Jl approximately 1. 0 or less. 

The experimental data for geometries C and D, which have large radius ratios are 
shown in figures 4(g) and (h) (p. 18), respectively. Noncrowbarred operation was ob- 
served for both with the glass baffle (shown by dashed lines in table I, p. 13) in place. 
Spectra indicated that both geometries had substantial quantities of electrode material in 
the exhaust. This erosion could affect the Jl chosen, especially in geometry D. Typi- 
cal plasma spectra are discussed in appendix C. No experimental velocity data were re- 
corded for geometry D since erosion seriously affected the data. 

Geometry D yielded substantially higher efficiencies than geometry C even though 
peak theoretical efficiencies predicted by the loss are almost the same (27. 8 and 28. 6 
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percent, respectively). Geometry D had a shorter barrel length to reduce wall losses, 
and the radius ratio was doubled compared with geometry C to show the effect of L* on 
performance while maintaining almost the same L*£ product. 

For geometry D, the maximum experimental efficiency is 22.2 percent, which is 
about 60 percent of the peak value predicted by the simplified model and about 7 5 percent 
of the peak value for the theoretical loss model. For 15 kilovolts, if Jl were increased 
to account for eroded electrode material in the exhaust, these experimental peaks might 
be near the predicted optimum values of Jl. 

The uncertainty in determining experimental velocity when probe traces are difficult 
to interpret is probably the reason for generally lower than predicted velocities especially 
for Jl greater than 1.0. In all but one case (geometry A, delayed internal crowbar), the 
velocities measured are significantly lower than the predicted values for Jl smaller than 
1. 0. These unexpectedly low velocities are probably the result of a compounding of ef- 
fects due to uncertainty in measuring velocity and Jl. 


CONCLUDING REMARKS 

Experimental calorimetric exhaust efficiencies and magnetic front velocities were 
compared with the prediction of a theory based on the snowplow model. Four gun geome- 
tries were studied for a range of initial masses of argon propellant and initial voltages. 
The guns were operated in several modes: noncrowbarred, crowbarred internally at the 
time of maximum current, crowbarred internally at a later, more nearly optimum, time, 
and crowbarred externally at maximum current time. 

Comparisons of the experimental results with the predictions of the theory indicate 
general qualitative agreement of efficiency and sheet velocity. The predicted kinetic en- 
ergy efficiencies and velocities, however, are consistently higher than measured values. 
These discrepancies might be attributed to an incorrect assumption for the mass distribu- 
tion because no measurements of distribution were made for this geometry. 

The predictions regarding the effect of crowbar were also in qualitative agreement 
with the experimental results. For one gun geometry, however, the peak efficiency for 
the crowbarred mode was larger than that of the noncrowbarred mode, contrary to the 
predictions of the theory. However, the presence of the glass baffle in the gun that oper- 
ated in the noncrowbarred mode might have significantly changed the shape of the mass 
distribution and thus obscured the comparison. 

The theoretical model presented is useful for the determination of relevant parame- 
ters and is adequate to predict qualitative trends. More accurate predictions probably 
require a more sophisticated model. 
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Because theory suggests the use of a sharply decreasing (slug) mass loading distribu- 
tion, future experimental efforts should be directed toward producing and measuring such 
distributions while attempting to suppress the instabilities that seem to be associated with 
these distributions. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 18, 1966. 
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APPENDIX A 


SYMBOLS 


C capacitance, F 

e electronic charge, 1.602x10"^®, C 

f nondimensional mass distribution 

function, M(x)/Mq 

g fraction of ions in excited state 

I current, A 

J nondimensional gun current, 

ig/k/ve^) 

j nondimensional heat conduction 

parameter (eq. (6)) 

K coefficient of thermal conductivity, 

J/(sec)(m)(°K) 
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k Boltzmann constant, 1.38x10 , 

J/°K 

L inductance, H 

Lq gun inductance from crowbar dis- 

charge location to initial dis- 
charge location, = L'jZq, H 

n 

L f inductance per unit length, 2x10” 

ln( r 0 / r j) , H/m 

& nondimensional inductance, 

(1/L*£)L 

i effective gun length, m 

M mass swept up by moving dis- 

charge, M = M(x), kg 

M q total injected mass per shot, kg 

Jl nondimensional mass loading pa- 

rameter, ^i/l'Qq^Mq 
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m g electron mass, 9. 107x10 kg 
m. ion mass, kg 

P nondimensional thermal ionization 

factor 

Q charge in capacitor, C 

Qq initial charge in capacitor, C 

R resistance, ohm 

0t nondimensional resistance, 

(C/L , £) 1 / 2 R 

r. radius of gun center electrode, m 
r Q radius of gun outer electrode, m 

S area, sq m 

T mean plasma temperature, °K 

t time, sec 

t mean lifetime of excited state, sec 

V capacitor voltage at time t, V 

Vq initial voltage supplied to the capac 
itor, V 

x distance of current sheet from ini- 
tial discharge position, m 

Ax thickness of current sheet, m 

z charge number of ions 

a nondimensional mass distribution 

parameter (eq. (7)) 

j3 nondimensional temperature, 

2k /L f C\ T 

m.V H ) 


24 



r ion diffusion current, particles/sec 

y nondimens ional charge, Q 

6 nondimensional diffusion coefficient, 
[l/2(r 0 - r,)]£ 

e Q permittivity of free space, 

8. 85xl0 -12 F/m 

7j nondimensional kinetic efficiency, 

l/2Jik 2 

A ratio of Debye length to Rutherford 

impact parameter 

4 nondimensional length, £~*x 

t nondimensional time, ( h'C£)~^^t 

4> e effective excitation potential per 

ion, J 

effective ionization potential per 
ion, J 

<f>j combined effective potential, 

+ g$ e , J 

$2 excitation loss rate parameter, 

$ 2 = g$ e A e , J/sec 

cp^ nondimensional effective potential, 

(Pi = $ 1 /(m.£/L'C) 

nondimensional excitation loss 

parameter, 

cp 2 = $ 2 VL’jfc/im.jf/L'C) 


£2 nondimensional convection cooling 

coefficient, 3/2 (m-l/L’C) j eVQ 

Subscripts: 

B bank 

C crowbar 

D from crowbar discharge location to 

initial discharge location 

G gun 

P plasma 

p parasitic 

tot total 

Superscripts: 

first derivative with respect to non- 
dimensional time 

second derivative with respect to 
nondimensional time 

* first derivative with respect to non- 

dimensional distance 
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APPENDIX B 


EXTERNAL CROWBAR SYSTEM 
by Fred F. Terdan 

"Closing” an ignitron switch when the voltage across it is only 100 volts or less in- 
troduces serious variation in closing time (jitter). This situation occurs when an ignitron 
is used as an external crowbar switch for a coaxial plasma gun. The gun should be crow- 
barred at the time of maximum current. This maximum current corresponds to near zero 
gun voltage, and at that instant, an ignitron crowbar switch will have excessive jitter. 

The desired approach was to use ignitron circuitry, but conventional ignitron circuitry 
will not switch controllably at zero voltage. An ignitron biasing technique was developed, 
and this technique allowed controllable switching to occur at zero voltage. 

The external crowbar systems described in references 23 to 25 were not practical for 
this experiment when studied in terms of the characteristics of the existing capacitor bank. 
The biased ignitron circuitry developed for this experiment was suitable in all character- 
istics with the exception that it was inefficient because the switch impedance was higher 
than the bank impedance. 

The schematic diagram for the external crowbar system is shown in figure 5. The 
complete unit consists of the crowbar switch (three pairs of parallel GL-7703 ignitrons), 
a trigger amplifier (incorporating three stages of amplification utilizing a 2D21 thyratron, 
a 5C22 thyratron, and a GL-7703 ignitron), a control system, and the necessary power 
supplies. Three back-to-back pairs of parallel ignitrons were chosen to provide neces- 
sary voltage holdoff capability, current carrying capacity, low inductance, and low volt- 
age switching ability. Capacitors of 1 microfarad placed across each set of three igni- 
trons were precharged to opposite polarities. Charging is accomplished by a remotely 
operated power supply just prior to the use of the crowbar switch. The charged capaci- 
tors provided some of the energy lost in making the crowbar switch conduct. These capac- 
itors and a water system that controlled ignitron temperature were necessary to reduce 
the switch jitter to less than 0. 2 microsecond. Figure 6 shows a view of the crowbar ig- 
nitrons, the capacitors in the crowbar system, and the gun header. The total inductance 
of the crowbar switch is 39 nanohenries, and the static crowbar switch resistance is ap- 

_3 

proximately 3. 5x10 ohm. This assembly holds off 25 kilovolts and can be activated 
(with 0. 2 jusec jitter time) for header voltages above a few hundred volts. 
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Figure 6. - External crowbar system. 
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APPENDIX C 


SPECTROSCOPIC OBSERVATIONS 

The discharges of two different gun geometries were examined with a 1^ -meter grat- 
ing spectrograph. The first-order dispersion of the instrument was 20 A per millimeter. 
Since the discharge is of such short duration, a series of 25 shots was superimposed to 
expose the spectrographic film sufficiently. Reflected light did not have sufficient inten- 
sity to record on the film. The spectra recorded are due to the time- integrated light 
from the discharges that were in the viewing area. Film limitations placed the lower 
viewing limit at 2400 A. 

The discharge in the annulus of geometry A was viewed radially (i. e. , normal to the 
gun axis) through a fused quartz window in the outer barrel. Figure 7(a) shows the spec- 
trograph obtained with geometry A operated internally crowbarred at 20 kilovolts and 
with argon propellant. (These spectra were identified and interpreted by Henry J. Hettel 
of Lewis.) The window was located 41. 5 centimeters from the breech end of the gun 
(6. 5 cm upstream of the gun muzzle). Argon n and argon HI lines are noted and the dis- 
charge is relatively free from contaminants. From these qualitative measurements, the 
effective ionization potential per ion for argon was estimated to be 30 electron volts. This 
value was used to calculate the nondimens ional effective potential <p-^ listed in table n 
(p. 19). The effective excitation potential was estimated as 40 electron volts per ion and 
was used to calculate the nondimensional excitation loss parameter ip^, also listed in 
table n. 

In contrast, the spectrum for a modified geometry D gun (radius ratio, 7. 5 and with- 
out the glass baffle), using argon is shown in figure 7(b). The device was also operated 
internally crowbarred at 20 kilovolts. The spectrograph viewed the exhaust in a part of 
the gun annulus as it left the gun muzzle. The line of sight was at a 45° angle to the gun 
axis, looking into the gun muzzle through a heat-resistant glass drift tube. Electrode 
material (copper and solder) and glass dominate the spectrum and no argon propellant 
lines are observed. The glass spectrum probably is caused by the excitation of the heat- 
resistant glass tube as the plasma exhaust collides with the tube wall downstream of the 
gun muzzle. The presence of copper lines and the lack of argon lines of any appreciable 
intensity indicates that the device is probably operating on metal- derived plasma or some 
combination thereof. 
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(a) Radial view near muzzle annulus of geometry A. 
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(b) Modified geometry D viewed at 45° angle to muzzle annulus. 

Figure 7. - Section of visible spectrum of coaxial plasma-gun discharge. 
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